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The LHC experiments have reached a milestone in our understanding of nature. The 
evidence for the observation of the Higgs spin-0-boson as a manifestation of a scalar field 
provides the so far missing corner stone for the standard model of particles and fields, "which 
describes the fundamental particles from which every visible part of the universe are made 
of and the forces between them" p~H3]. However, the standard model fails to explain the 
other, non-visible but gravitationally active part of the universe, contributing to about 95% 
of the dark energy/mass density of the universe. Its nature is unknown but the confirmation 
of a scalar Higgs is giving a boost to scalar-field-theories in an attempt to explain the other, 
non-visible part of the universe. So far gravity experiments and observations performed at 
different distances find no deviation from Newton's gravity law and also fail to find a fifth 
force contributed to dark energy or dark matter 0]. Therefore dark energy must possess a 
screening mechanism which suppresses the scalar-mediated fifth force. Our line of attack is 
a novel gravity experiment with neutrons based on a quantum interference technique. In 
particular the spectroscopic measurement of quantum states on resonances with an external 
coupling makes this search in the gravity potential of the Earth a powerful way to search for 
dark matter and dark energy contributions in the universe. Quantum states in the gravity 
potential are intimately related to other scalar field or spin-0-bosons if they exist. If the 
reason is that some undiscovered particle interact with a neutron, this should result in a 
measurable energy shift of quantum states in the gravity potential of the Earth, because for 
neutrons the screening effect is absent [5]. Here we use Gravity Resonance Spectroscopy [6] 
to measure the energy splitting at the highest level of precision reported so far, providing a 
constraint on any possible new interaction. We obtain a sensitivity of 10 -14 eV. A fit to the 
data sets an upper limit on any fifth force, in particular on parameter /3 < 2 x 10 9 at n = 3 
for the scalar chameleon field, which is improved by a factor of 100 compared to the upper 
bound 5 derived from our previous experiment [7J [H| and five orders of magnitude better 



than the upper bound, j3 < 10 14 , from precision tests of atomic spectra [5]. The pseudoscalar 
axion coupling is constrained to g s g p /hc < 3 x 10~ 16 at 20 pm, which is an improvement 
by a factor of 30. These results, consistent with zero, indicate that gravity is understood at 
this improved level of precision. 
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The last fourteen years of experiments and observation [9] have produced an entirely new map 
of the universe. Cosmological parameters describe the properties and the global dynamics of the 
universe. Observational cosmology has determined these parameters to one or two significant figure 
accuracy. Of great interest is now how the energy-matter budget is built up from its constituents: 
e.g. baryons, photons, neutrinos, dark matter and dark energy [10J. The known particles of the 
Standard Model |11H13| account for only 5%, whereas the majority consists of unknown dark energy 
and dark matter. 

Firstly, dark energy - interpreted as negative pressure - is causing the universe to accelerate. A 
particularly attractive possibility is that a scalar field, usually called quintessence, is responsible for 
inflating the universe |14H17j . A competing possibility is that the acceleration is due to a modifi- 
cation of gravity, i.e., the left-hand side of Einstein's equation rather than the right. Observations, 
which can distinguish these two possibilities are desirable, since measurements of expansion kine- 
matics alone are not able to. Remarkably, the best limits for scalar fields (chameleons) come from 
our experiment with ultra-cold neutrons, as we will show. 

Secondly, distinct from dark energy is the concept of dark matter observed indirectly, through 
the influence of its gravity on surrounding ordinary matter. It has evaded more conventional 
observation using optical, radio or other telescopes, so it evidently only has a feeble coupling to 
photons. We do not yet know precisely what this dark matter is, but the leading theoretical 
candidates (axions and weakly interacting massive particles, WIMPs) are massive particles that 
do not acquire their mass by interacting with the electroweak Higgs condensate. One eye-catching 
concept is the idea of detecting very light bosons through the macroscopic forces which they 
mediate. These particles may be responsible for the dark matter in the universe or part of it. 
They must be very light with a mass of about 10 ueV to 1 eV to be compatible with astronomical 
observation [18]. This corresponds to a range of 0.2 um and 2 cm they can mediate, leading to a 
deviation from Newton's law at short distances, exactly in the range of this experiment. 

Thirdly, there are indications that all known basic forces, the weak force, the strong force, 
electromagnetism and gravity, might be unified at high energies. One approach are string-theories, 
which are naturally defined in ten or eleven space-time dimensions. The hope is that a future theory 
of quantum gravity can be constructed in such a way that all forces, including gravity, could be 
combined within it. Such a theory cannot be constructed in a consistent way if space-time is 
limited to four dimensions, so additional spatial dimensions are needed to accommodate such a 
theory in a reasonable way. Extra-dimensions have not been observed in everyday life, because 
these extra-dimensions might be "curled-up" with a possibly very small radius of curvature. This 
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should lead to deviations from Newton's gravitational law at very small distances. The focus of 
these has shifted over the years and it has been realized, that the string scale, which controls the 
strength of gravity as well as the strength of other forces, is not necessarily related to the Planck 
length of 10 -35 m, instead it can be considered as a free parameter and we might live with large 
extra dimensions |19j . again in the range of our experiments |20j. 

The idea of finding a common ground to the problems related to string theories, dark energy and 
dark matter is aesthetically appealing to us. That these problems might have a common solution 
can be seen that these approaches are strongly interrelated, not only by a unified physics aim, like 
the role of gravity in modern gauge theory for cosmology, but also by the common methodical 
approach addressed by this experiment. 

All these ideas triggered gravity experiments of different kinds, which in the past ten years have 
validated Newton's gravitational law down to about 20 um [4, 20-28J. Searches for short-ranged 
non-Newtonian gravity so far have been performed using neutrons, atoms, torsion pendulums, and 
other macroscopic objects. To date, neutron scattering experiments provide most stringent limits 
at very short distances around 1-10 nm \29\ [30J. A recent review on "the neutron and its role in 
cosmology and particle physics" can be found in Ref. |31j . Currently, measurements with atomic 
force microscopes determine the limits on non-Newtonian gravity below 10 um. 

Such experiments clearly show the basic problem in searching for new physics at small distances: 
that the size of the objects under study needs to be reduced, which is accompanied by a reduction 
in signal intensity. At the same time, the electrostatic background increases. The gravitational 
interactions are studied in the presence of large van-der-Waals and Casimir forces, which depend 
strongly on the geometry of the experiment, and the theoretical treatment of the Casimir and 
other polarizability effects is a difficult task. In practice, the experimental data are subject to 
corrections, which can be orders of magnitude larger than the effects actually searched for [3]. 

Our solution to this problem is the development of a purely quantum mechanical system for 
gravity measurements, where polarizability effects for all practical reasons are absent, i.e. the 
excitation of bound quantum states of a neutron in the gravity potential of the Earth, and its 
interaction with a macroscopic system, here a mirror for neutron reflection. This so-called Gravity 
Resonance Spectroscopy technique (GRS) does not rely on electromagnetic fields or a coupling to 
an electromagnetic potential [32] . 

To improve on this resonance techniques we have developed a novel frequency reference system 
for gravity measurements, see Fig. [T] This system has two great advantages: Firstly, the reference 
system is only based on natural constants: the mass of the neutron m, the Planck constant h, 
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FIG. 1. Left: Frequency reference for gravity experiments based on energy eigenvalues of a neutron in the 
gravity potential of the Earth with confining mirrors at the bottom and at the top separated by 30.1 pm. Also 
shown is the neutron density distributions for energy levels |1) to |5). The Gravity Resonance Spectroscopy 
technique allows to drive transitions between these states. Quantum transitions |1) ■<-> |2), |1) ■<-> |3), 
|2) O |3), and |2) -H- |4) have been observed. 

Right: Resonant transition frequencies in Hz at width I = 30.1 pm between mirrors, see Methods. 

and the acceleration of the Earth g. The constants m and h are known within and g can 

be measured with a corner cube within Ag/g = 2 x 10~ 9 . Tidal changes affect g within 10~ 7 . 
Such a notional level of accuracy for the energy spacing in the gravity potential of the Earth is a 
challenge for our experiment to be probed by GRS. It corresponds to an energy change of around 
AE = 10 -20 eV. On this level of accuracy, we want to apply GRS to the measurements of the 
inverse-square law of the Newton gravity at length scales of the order of a micron, where it is a solid 
testing-ground for theoretical schemes of interactions, based on supersymmetric theories of strings 
and D-branes, predicting possible deviations from Newton's Inverse Square Law. The lowest energy 
eigenvalues E n , (n > 1) are on the pico-eV energy-scale. In this experiment these energy eigenvalues 
are tuned by a second mirror on top of the bottom mirror at a height of 30.1 pm as described in 
Ref. [6]. This changes the resonance frequency between \p) and \q) due to the additional mirror 
potential, which shifts the energy of states n > 1, but leaves state |1) unchanged. The additional 
potential provides an effective increase in sensitivity by a factor of two in a search for hypothetical 
axion or chameleon induced phase shifts or for other fifth forces. Secondly, we have now on hand 
precision measurements of resonant transitions between states |1) «-)■ |2), |1) o |3), |2) «-)■ |3), and 
1 2) «-» 1 4) in the frequency regime between 70 Hz and 800 Hz. This quantum system represents a 



6 



1.2 




0.4 b — ' — ' — ' — 1 — ' — ■ — ■ — 1 — ■ — ■ — ' — 1 — ' — ' — ' — 1 — — r I , , , , I , , , , I , , , , I , , , , I , , , , I 

200 400 600 800 1 2 3 4 5 6 

Frequency [Hz] Vibration amplitude [mm/s] 



FIG. 2. Gravity Resonance Spectroscopy: 

Left: Transmission as a function of frequency. Transitions between |1) <-> |2) , |1) -f-> |3), |2) •<->• |3), and 
1 2) O 1 4) have been observed. We find transmission frequencies z/12 = (258 ± 2) Hz, V23 = 280^2 Hz, 
v 13 = (539 ± 3) Hz, v 2i = 6791^ Hz. 

Right: Damped Rabi-oscillations and Transmission as a function of amplitude strength in resonance and at 
two other positions in the rising slope, see arrows to the left. The resonant frequencies have been determined 
to up to 5 x 10~ 3 corresponding to an energy sensitivity of E = 10 -14 eV. 



self-consistent frequency system, where systematic effects are virtually absent. It can be compared 
- however with much less precision so far - with a frequency reference in laser or microwave 
spectroscopy, e.g. atomic clocks, nuclear magnetic resonance techniques or Mofibauer spectroscopy, 
where the observables have generally been restricted to electromagnetic interactions. Our current 
sensitivity is E = 10~ 14 eV, see Table|and Fig. g 

Fig. [4] shows the apparatus consisting of a mirror underneath and a mirror on top at distance I 
with rough surface. The neutrons are taken from the ultra-cold neutron installation PF2 at Institut 
Laue-Langevin (ILL). The experiment itself is mounted on a polished plane granite table with a 
passive anti-vibration tables underneath. The granite is leveled to an accuracy of better than 
10 urad and stabilized with a precision better than 1 urad. A mu-metal shield suppresses the 
coupling of residual fluctuations of the magnetic field to the magnetic moment of the neutron 
sufficiently. Gravity Resonance Spectroscopy is a Rabi-type spectroscopy method and consists of 
a state selector, where the initial state \p) is prepared. Transitions into a second state \q) are 
induced by applying a so-called 7r-pulse. A detector at the end measures these neutrons, which 
were accepted by the state selector. On resonance a sharp decrease in count rate is found. The 
frequency width of the resonance is given by the time the neutron spends in the oscillating system. 

In general, the oscillator frequency at resonance for a transition between states with energies 
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E p and E q is given by 



1 Eg- E p 

2vr h ' 1 J 



The linear gravity potential leads to discrete energy eigenstates of a bouncing quantum particle 
above a horizontal mirror. Neutrons are trapped by gravity, but the motion is unrestricted parallel 
to the mirror. So far the systematic effects of the top mirror on the resonant transitions are smaller 
than our statistical error, see Methods. A measured energy level splitting can be compared with 
the prediction in the gravity potential, allowing a test of Newton's law at short distances. The 
lowest energy eigenvalues E n read 1.41 peV, 2.50 peV, 3.66 peV, 5.22 peV, and 7.25 peV. Another 
GRS measurement, based on oscillating magnetic gradient fields, is in progress |33| . 

Our measurement is derived from 116 neutron transmission measurements taken in 2010 and 
2011, together with many subsidiary measurements used to search for systematic errors. The 
neutron transmission in counts per second as a function of frequency and amplitude is registered. 
Included are 17 measurements, where the vibration amplitude is zero. After background subtraction 
of (2.18 ± 0.08) x 10 -3 s _1 , the following 10 parameters fit the 116 data points, using Eq. 
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Tab.[Ij Firstly, three parameters 172 ; 1^23, and ^24 correspond to the transition frequencies. Energy 
conservation requires Z72 + ^23 = ^13 • Secondly, three parameters 72, 73, and 74 describe the 
damping of the probability amplitude of the states due to the rough, upper surface of the neutron 
mirror with respect to the damping of the ground state 1). Thirdly, the initial populations C2, C3 
and C4 for state 2) to 4) with respect to state 1) and an overall normalization. In a first analysis, 
resonant transitions 172, ^23, and z/24 are determined, Z73 is obtained from energy conservation. A 
Xmin = H4.7 at 116 data points and 10 parameters is found, which corresponds to a probability of 
26.6%. Resonant frequencies are determined with a precision of 8 x 10~ 3 for transition 1 •<-)■ 2, 1% 
for transition 2 f> 3, 5 x 10~ 3 for transition lf>3, and 2% for transition 2 f> 4, see Tab.|lj In fact, 
the resonant transitions are given by the acceleration of the Earth and the slit height I between the 
two mirrors. Due to the reduced number of parameters, the \ 2 is worse, but with g = 9.81 m/s 2 , 
which is the value for Grenoble, the probability is still 16.7%. 

Fig. [2j left, shows the measured relative transmission. Clearly visible is a decrease of transmis- 
sion as a function of frequency. Such dips correspond to relevant transition frequencies Z72, Z73, 
f23, and z^24- At a frequency v = 280 Hz, the absolute resonance minimum is found. Scanning 
the vibration amplitude generates a damped Rabi-oscillation curve as predicted, shown in Fig. [2j 
right. These measurements altogether can be displayed in one single data point corresponding to 
one specific frequency (black arrow in Fig. [2j left), by normalizing the measured count rate to the 
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TABLE I. Best parameter set at Xmin anc ^ error at la standard deviation (Xmin + -0 m ^ dimensional 
hyperspace. 

average rate without vibration multiplied by the best fit for a chosen vibration strength divided by 
the best fit at that point. The instruction offers an elegant way in displaying all 116 measurements 
in one Figure. [2j left, irrespective of vibration strength. For these relative transmission points a 
binning of 20 Hz is used. The chosen vibration strength is 2.03 mm/s, which is the average of all 
vibration measurements. 

In this section we perform a search for so called chameleon fields. A chameleon field has been 
suggested to drive the current phase of cosmic acceleration for a large class of scalar potentials. The 
properties depend on the density of a matter to which it is immersed. Because of this sensitivity 
on the environment, the field was called chameleon, and it can couple directly to baryons with 
gravitational strength on Earth but it would be essentially massless on solar system scales |34H37] . 
The coupling constant is usually referred to as (3. Such a chameleon coupling to our neutron would 
provide a frequency shift of resonant transitions, which is proportional to (3 and can be measured by 
comparing the transition frequency v vq with its theoretical expectation fpg eo , when no chameleon 
field is present. Our method directly tests the chameleon matter interaction and does not rely on 
the existence of a chameleon-photon-interaction as other experiments do |3 81 - 133] . 

The calculation of the contribution of a chameleon field to the transition frequencies of the 
quantum gravitational states of ultra-cold neutrons is related to the calculation of the matrix 
elements {p\(ft(z)), where \p) is a low-lying quantum state of ultra-cold neutrons in the gravitational 
field of the Earth. Due to the Airy functions describing the gravitational states of ultra-cold 
neutrons, the main contribution to matrix elements {p\z\p) comes from the region around (d/2 + 
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FIG. 3. Exclusion plots with 68%, 90%, and 95% C.L.. Left, chameleon field: limits are shown as contours 
for parameter f} vs. parameter n, which is improved by a factor of 100 [5] compared to the upper bound 
derived from our previous experiment 018] and 5 orders of magnitude better than the upper bound, j3 <10 14 , 
from precision tests of atomic spectra [5]. Right, axion field: limits are shown as contours for parameter 
g s g p /hc vs. range A. A limit for an attractive coupling is shown in a solid line, whereas a repulsive force 
limit is indicated by the dashed line. 

z) Iq = (2m 2 g/h 2 )^ 1 / 3 = 5.9 um, which is a natural scale for neutron quantum states. In our case 
with two mirrors at distance d, the chameleon field potential is given by [44J: 

2 

(d\ m ( 4z 2 \ "+ 2 

2 

(n -L- 9 \ n+2 

^ Ad ) > ( 3 ) 

where Mpi is the Planck mass, the scale A is chosen to be equal A = 2.4(1) x 10~ 12 GeV. Quanta of 
a chameleon field are massive with mass defined by the second derivative of the potential V((j>). The 
coordinate system has been shifted by d/2 due to symmetry reasons. In a fit to the data, the Earth's 
acceleration has been fixed at g = 9.81 m/s 2 , the other parameters have been varied. The extracted 
confidence intervals for limits on the parameters (5 and n can be found in Fig. [3} The experiment is 
most sensitive at n=3 and excludes a coupling to chameleon fields with /3 > 2.1 x 10 9 (90% C.L.). 
This is a factor of 100 more precise than the previous upper limit [5]. 

Second, we perform a search for spin-dependent forces at short distances. Resonant transitions 
are for this purpose measured with polarized neutrons. The experiment is therefore slightly mod- 
ified: a homogeneous magnetic guiding field of 50 uT preserves the neutron spin throughout the 
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experiment. Neutron spin polarization is analysed by a soft ferromagnetic foil, which was coated 
on the detector foil for that purpose. A strong magnetic field at the position of the detectors is 
oriented parallel or antiparallel to the direction of gravity. By altering the direction of this field 
the analysing power with respect to the spin state direction and gravity direction can be changed 
using the hysteresis of the iron. A hypothetical spin-dependent force would change the transition 
frequencies. Such a frequency shift is measured by reversing the direction of the applied guiding 
field and the detector field and measuring the difference in the two signals. As a check, the signal is 
measured at two transition frequencies. This difference is positive on one side of the transmission 
curve and negative on the other. The most sensitive places are those with large slopes. We do not 
observe a significant frequency shift. Limits for spin-dependent forces can be easily interpreted as 
bounds of the strength of the matter coupling of axions. Axion interactions with a range within 
0.2 um < m a < 2 cm (corresponding to axion masses (10~ 5 < m a < l) eV), the " axion window", 
are still allowed by the otherwise stringent constraints posed by cosmological data. An axion would 
feel a CP-violating spin-dependent interaction in the presence of matter given by [45] . 

Here, a denotes the neutron spin and ft is a unit vector related to the geometry of the macroscopic 
matter configuration. 

A fit of strength g s g p /hc and range A together with the other independent parameters, this 
leads at 90% C.L. to an upper limit on the axion interaction strength as shown in Fig.|3j right. For 
example, at A = 20 um in the astrophysical axion window, an attractive coupling strength g s g p /hc < 
3 x 10~ 16 is excluded. This is the most precise upper limit from a direct search for attractive and 
repulsive g s g p / he-coupling at this length scale A. It is a factor of 30 more precise than a limit 
that has been derived [36] from our previous experiment with ultra-cold neutrons [3 |S]. Other 
experiments have searched for spin-mass coupling at larger A and extrapolated over several orders 
of magnitude. A recent review |31j is covering this topic. The most restrictive limit on the product 
gsgp/hc has been derived by combining the existing laboratory limit on the scalar interaction g s 
with stellar energy- loss limits on the pseudoscalar coupling g p , which is more restrictive than 
laboratory searches |47] . 

Our method profits from small systematic effects in such systems, mainly owing to the fact that, 
in contrast to atoms, the electric polarizability of neutrons is extreme low. Our error is dominated 
by the statistical uncertainty. Our experiment leads the way in the application of GRS to probe 
for new particle physics and to search for Non-Newtonian gravity with high precision. 
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FIG. 4. Left: a schematic diagram of the GRS setup. 1: bottom neutron mirror; 2: top rough mirror; 
3: reference mirror; 4: micro-positioning table; 5: piezo-control; 6:adjustment; 7: fine pitch thread; 8: 
spacer. Right: neutron density distribution above the mirror in front of the detector, measured by a track 
detector 48 51j. A fit to the data (red) is compatible with a population of 67% for state |1) (grey) and 33% 
for state |2) (green). 



METHODS 



The setup is shown in Fig. |4| Gravity Resonance Spectroscopy is a Rabi-type spectroscopy 
method and consists of a state selector, where the initial state \p) is prepared. By applying a 
so-called 7r-pulse, transitions into a second state \q) are induced. A detector at the end measures 
these neutrons, which were accepted by the state selector. On resonance a sharp decrease in count 
rate is found. The frequency width of the resonance is given by the time the neutron spends in the 
oscillating system. In our case the resonant transitions |1) <-> |2) and |2) -H- |3) are not separated; 
there exists a frequency band, where transitions |1) |2) |3) can simultaneously be driven, but 
they are in frequency space well separated from |1) f-> |3) transitions. The setup is consisting of a 
polished mirror on the bottom and a rough scatterer on top. In addition a third mirror serves as 
reference. Vibrations are controlled with a three-beam laser interferometer. 

A neutron entering this system is in a superposition of different eigenstates (p n . The wavefunc- 
tion as a solution of the Schrodinger equation between the mirrors with distance I, acceleration g 
and height z is given by 

4, (z, l,g) = J2 cne-^e-^^n (z, I, g) . (5) 

n 

c n is the state population of state n at the entrance of the device, the damping term e _7n//2 ' takes 
the damping of state n during the time t of flight into account, where the factor 7„ depends on 
the roughness of the upper mirror, and e~ l ^ n is a phase factor. When the vibrations are turned 
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FIG. 5. Top: FFT of a typical frequency measurement with a sampling rate of 2604 s _1 . The frequency 
is 280 Hz at a vibration amplitude of 2.12 mm/s. Resonant frequencies of the mirror setup at 120 Hz and 
of the turbomolecular pump at 833 Hz are visible. Left: Inverse FFT and a fit to vibration scan. Right: 
Residuals for amplitude measurement. 

off, the population in front of the detector is 67% for state |1), and 33% for state |2), state 
1 3) and higher states are zero. The state selector is a mirror with a rough surface of roughness 
3 um, and suppresses the population of higher levels by scattering and absorbing them. At the 
beginning of the experiment, the state population has been extracted from \ij)(zj, g)\ 2 , the spatial 
neutron distribution above the mirror i.e. measured with a track detector |49H51j . see Fig.|4j right. 
\ip(z, l,g)\ 2 is an incoherent superposition of state |1) and |2) since the phase is not defined. For the 
spectroscopy measurements described in this paper, this spatial resolution detector is replaced by 
a counting detector. In that detector, the entrance window is made out of aluminium coated with 
10 B, converting the neutron to 7 Li and a with an efficiency of about 92%. The background rate is 
Nq = (2.18 ± 0.08) x 10~ 3 s _1 . Background measurements are done every 200 s and in addition six 
longtime measurements are performed. The signal is registered with a usual preamplifier, main- 
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amplifier and multi-channel analyzer chain. In addition, a SADC-analyzer measures the identical 
neutron signal, allowing a post procession of the data due to energy and time registration as a 
function of time. For spin dependent searches, the detector foil is replaced by a foil coated with 
soft iron on the other side serving as an analyser for the direction of the spin of the neutron. The 
degree of the foil polarization has been measured separately to be ca. 93%. During all measurement, 
a beam monitor registered the neutron flux. The horizontal neutron velocity is (6 < v x < 10) m/s, 
which was adjusted with baffles in front of the setup. 

For the Rabi spectroscopy measurement, the observable is the measured transmission as a 
function of the modulation frequency and amplitude. This is done using a vibrating mirror i.e. 
a modulation of the mirror's vertical position. For this purpose, piezo-elements are attached 
underneath, see Fig [4} They induce a rapid modulation of the surface height with amplitude A. 
In accordance with transmission measurements of Fig. [2j resonant transitions are possible between, 
|1) «-)■ |2), |2) -H- |3) ,|2) «-)■ |4), and |1) «-)■ |3). The resonance width of about 53 Hz causes 
simultaneous transitions between |1) <-> |2) and |2) |3), they "overlap". As a consequence the 
transitions are treated as three-level systems and the time evolution of the states are given by the 
following differential equation: 
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(8) 

Here, for two states \p) and \q) with relative phase 4> qp , S qp is given by 

S qP = \e- i ^e- i ^ t d{q\^ z \p). (9) 

The frequency detuning from resonance reads 5 qp = 2tt (v — v qv ). The vibration amplitude in units 
of a velocity is given by u = 2irvak- 

The first equation examines transitions |1) |2) and |2) f-)- |3) for states |1), |2), |3), whereas 
the second equation system examines |1) |3) and |2) o |4) for states |1), |2), |3), and |4), which 
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is simpler because these transitions do not overlap. 

Vibration control is done in a twofold way. Firstly, a four-channel acceleration sensor-system 
measures external vibrations at different positions for surveillance. Secondly, a three-beam laser 
interferometer is in use. Laser beam LI serves as a reference beam, which measures vibrations on 
a separate mirror placed directly on the granite table. Laser L2 scans the surface of the vibrating 
neutron mirror to yield a spatial- and time-dependent map. The total deviation of the vibration 
amplitude from its mean was usually less than 20%. L3 scans for vibrations at a fixed position on 
the neutron mirror. From this data, a vibration analysis to extract the vibration amplitude a,k and 
relative phase at different positions of the mirror surface is performed using the expression 



where Nq is a constant offset, amplitude a/-, for frequency v\. and Phase <pk- The number of fit 
parameters are 3K + 1, where k is the number of sine- modes. 

In all measurements, a linear gradient of the vibration amplitude on the surface was found, 
which has no influence on the Gravity Resonance Spectroscopy method to first order. 

A typical frequency measurement is shown in Fig[5j top. After fast fourier transformation (FFT) 
a gauss filter is applied and inverse FFT in order to keep computing time small for adjustment 
to the sine- function, see Fig. [5j left and right. A two-dimensional map of the vibrating mirror for 
amplitude and frequency is shown in Fig. [6} 
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FIG. 6. Two-dimensional map of the vibrating mirror for amplitude (left) and phase (right). 
A grid of 12 x 10 (blue) data points with each 2 13 data samples taken at a sampling rate of 2604 s _1 results 
in a frequency resolution of 10 . The left figure shows, that the vibration amplitude is homogeneous on 
the 10%-level. Its relative phase is close to zero, as the plot on the right indicates. The red labels indicate 
the laser beams L2 and Li. 
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